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CONSISTENT  PI  ANALYSIS  OF  AQUEOUS  URANIUM-235 
CRITICAL  ASSEMBLIES 
By  Daniel  Fieno 


SUMMARY 

The  lethargy-dependent  equations  of  the  consistent  PI  approximation 
to  the  Boltzmann  transport  equation  for  slowing  down  neutrons  have  been 
used  as  the  basis  of  an  IBM  704  computer  program.  Some  of  the  effects 
included  are  (l)  linearly  anisotropic  center  of  mass  elastic  scattering, 
(2)  heavy  element  inelastic  scattering  based  on  the  evaporation  model  of 
the  nucleus,  and  (3)  optional  variation  of  the  buckling  with  lethargy. 

The  microscopic  cross-section  data  developed  for  this  program  covered 
473  lethargy  points  from  lethargy  u  =  0  (10  Mev)  to  u  =  19.8  (0.025  ev). 

The  value  of  the  fission  neutron  age  in  water  calculated  here  is 
26.5  square  centimeters;  this  value  is  to  be  compared  with  the  recent 
experimental  value  given  as  27.86  square  centimeters.  The  Fourier  trans¬ 
form  of  the  slowing-down  kernel  for  water  to  indium  resonance  energy 
calculated  here  compared  well  with  the  Fourier  transform  of  the  kernel 
for  water  as  measured  by  Hill,  Roberts,  and  Fitch. 

This  method  of  calculation  has  been  applied  to  uranyl  fluoride  - 
water  solution  critical  assemblies.  Theoretical  results  established 
both  unreflected  and  fully  reflected  critical  assemblies  have  been  com¬ 
pared  with  available  experimental  data.  The  theoretical  buckling  curve 
derived  as  a  function  of  the  hydrogen  to  uranium-235  atom  concentration 
for  an  energy- independent  extrapolation  distance  was  successful  in  pre¬ 
dicting  the  critical  heights  of  various  unreflected  cylindrical  assem¬ 
blies.  The  critical  dimensions  of  fully  water-reflected  cylindrical 
assemblies  were  reasonably  well  predicted  using  the  theoretical  buckling 
curve  and  reflector  savings  for  equivalent  spherical  assemblies. 


INTRODUCTION 

The  NASA  Lewis  Research  Center  has  for  some  time  operated  a  Zero 
Power  Reactor.  This  reactor  consists  of  uranyl  fluoride  salt  (UOcFg) 


2 


dissolved  In  water  and  contained  in  a  cylindrical  aluminum  tank.  The 
parameters  of  the  system  are  the  diameter  of  the  cylindrical  tank,  use 
of  a  side  reflector,  and  concentration  of  uranyl  fluoride  salt  in  water. 
Criticality  is  achieved  hy  varying  the  height  of  the  solution  in  the 
tank.  No  control  rods  are  inserted  into  the  solution  while  the  reactor 
is  critical.  Thus,  good  experimental  conditions  exist  for  obtaining 
criticality  data  and  performing  reactivity  measurements  hy  measuring 
changes  in  the  solution  height  due  to  insertions  in  the  reactor.  Descrip¬ 
tion  of  this  reactor  and  its  various  characteristics  are  found  in  the 
final  hazards  report  (ref.  l)  prepared  for  the  Atomic  Energy  Commission. 

The  criticality  data  obtained  from  the  NASA  Zero  Power  Reactor  as 
well  as  that  obtained  from  extensive  experimental  studies  performed  at 
the  Oak  Ridge  National  Laboratory  (refs.  2,  5,  and  4)  can  be  analyzed 
by  various  theoretical  models.  Callihan  (ref.  5)  describes  the  use  of 
Fermi  age-diffusion  and  Goertzel-Selengut  analytical  models  as  well  as 
a  method  due  to  Greuling  for  predicting  critical  reactor  diameters  of 
various  Oak  Ridge  experiments.  A  more  recent  study  of  the  Oak  Ridge 
experimental  data  by  Gwin  and  others  (ref.  6)  makes  use  of  the  empirical 
kernel  method,  which  is  based  on  the  measured  slowing-down  distribution 
in  water  rather  than  on  any  theoretical  model  of  the  neutron  slowing - 
down  process.  In  addition,  the  use  of  the  Los  Alamos  DSN  transport  code 
for  obtaining  minimum  critical  dimensions  for  water  solution  is  described 
by  Mills  (ref.  7). 

The  objective  of  the  present  study  is  to  provide  a  theoretical  model 
based  or.  the  numerical  solution  of  the  lethargy-dependent  consistent  PI 
approximation  to  the  Boltzmann  transport  equation  for  slowing  down 
neutrons.  This  model  is  coded  for  the  IBM  704  electronic  computer  and 
employs  microscopic  nuclear  data  (refs.  8,  9,  and  10).  This  theoretical 
model  is  tested  by  comparison  with  the  experimental  age  and  Fourier 
transform  of  the  slowing-down  kernel  of  fission  neutrons  in  water.  This 
model  is  further  employed  in  estimating  the  absolute  reactivity  of  un¬ 
reflected  U02F2-water  systems  and  is  used  as  a  basis  for  evaluating  con¬ 
sistent  three-group  parameters  for  the  interpretation  of  water-reflected 
UOgFg-water  systems. 


The  next  section  gives  the  derivation  of  the  general  lethargy- 
dependent  consistent  PI  equations  used  as  the  basis  for  an  IBM  704  code. 
Next  follows  a  section  giving  the  analytical  results  based  on  this  model 
and  comparison  with  experimental  data  obtained  with  the  NASA  Zero  Power 
Reactor  as  well  as  data  obtained  at  Oak  Ridge. 


THEORY 

Slab  Geometry  Spherical  Harmonics  Form  of  Transport  Equation 

Let  n(r, v;dr  dv  be  the  probable  number  of  neutrons  in  the  space 
volume  element  dr  about  position  r  and  in  the  velocity  element  dv 
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about  velocity  v  per  unit  time.  Then,  the  Boltzmann  transport  equa¬ 
tion  for  sieving  down  neutrons  at  steady-state,  which  equates  the  losses 
per  unit  time  of  neutrons,  due  to  transport  out  of  and  collision  re¬ 
moval  from  the  phase  volume  element  d?  d^,  to  the  gain  per  unit  time 
of  neutrons,  due  to  elastic  and  inelastic  in-scattering  and  fissioning 
of  the  fissile  nuclei  in  the  phase  volume  element  d?  dv,  has  the  form 

Q  •  Vty(r,v)  +  Et(v)cp(r,v)  =  J  dv'2es(v'  )cp(r,  v'  )P(^’  -*•  v) 


+  y*  dvr£in(v‘ )cp(r,  vf  )h(v*  -*  v)  +  S(?,v)  (l) 

where  the  neutron  flux  is  defined  as  cp(r,  v)  =  vn(r, v)  with  v  being 
the  neutron  speed)  S(r,  v)  is  the  source  of  neutrons  due  to  fiBsion; 

St(v),  SeB(v),  £in(v)  are  the  macroscopic  total,  elastic  scattering, 

and  inelastic  scattering  cross  sections,  respectively,  at  speed  v;  ft 
is  a  unit  vector  along  the  direction  of  motion  of  the  neutron;  P£v*^-*  v) 
is  the  probability  that  a  neutron  in  the  phase  volume  element  dr  dv' 
will,  after  an  elastic  scattering  collision,  be  in  the  phase  volume 
dr  dv.  Similarly,  H(v'  -*  v)  is  the  probability  that  a  neutron  in  the 
phase  volume  element  d?  dv'  will,  after  an  inelastic  scattering  col¬ 
lision,  be  in  the  phase  volume  d?  17. 

Now  transform  equation  (l)  from  (?, v)  apace  to  (?,  v,  ft)  space  and 
finally  from  (tv  y.  o )  space  to  (r,u,  ft)  space  where  u  is  the  neutron 
lethargy  defined  by  the  relation  u  «  ln(Eg/E)  with  Eq  usually  taken 

as  107  ev.  The  result  of  these  transformations  and  writing  an  explicit 
form  for  the  fission  source  is 

3*  ^p(?, u,ft)  +  £t(u)9(r, u,tf)  =  /du'Ies(  u’)/  dll'  cp(r,  u;ff'  )P(u-u;  uL) 

+  J  duTiin(u')  y*  dB'cpCrjUjfl' ) 

+  (2) 

where  v(u)  is  the  number  of  neutrons  released  per  fission  occurring  at 
lethargy  u,  Ef(u)  is  the  macroscopic  fission  cross  section,  f(u)  is  the 

fraction  of  the  fission  neutrons  born  in  the  lethargy  element  du  about 
u  (the  integral  of  f(u)  over  all  lethargies  is  normalized  to  unity), 

Keff  is  the  static  criticality  factor,  and  dft  is  the  element  of  solid 

angle  about  the  direction  ft.  The  inelastic  scattering  kernel  is  assumed 
to  depend  only  on  the  initial  lethargy  u'  and  the  final  lethargy  u 
of  an  inelastically  scattered  neutron  and  to  be  isotropic  in  the  labora¬ 
tory  coordinate  system.  The  elastic  scattering  kernel  is  assumed  to 
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depend  only  on  the  difference  between  the  final  and  initial  lethargies 
and  on  the  cosine  of  the  angle  q^  between  the  initial  and  final  di¬ 
rections  of  motion  of  a  neutron  undergoing  an  elastic  scattering  col¬ 
lision.  Davison  (ch.  II,  ref.  ll)  gives  additional  details  concerning 
equation  (2). 

If  the  neutron  distribution  varies  only  along  the  Z-direction  and 
is  azimuthally  symmetric  about  that  axis,  if  equation  (2)  is  Integrated 
over  the  azimuthal  angle,  and  if  the  flux  cp  (=  cp(Z,  u,  q)  (the  quantity 
q  is  the  cosine  of  the  angle  between  the  Z-direction  and  the  direction 
of  neutron  motion  fi)  and  the  elastic  scattering  kernel  are  expanded 
into  a  series  of  Legendre  polynomials,  that  is, 

00 

<P(Z,u,q)  =  2  - % —  cPj(^u)P2(q)  (3a) 


?(u-u,*qI)5  a  ^  I"-1-  P^u-U*  )Pz(uL)  (3b) 

2=0 


and  these  expansions  are  substituted  into  the  equation  for  cp(Z,u,  q), 
which  is  then  multiplied  by  a  Legendre  polynomial  Pj^(q)  and  integrated 
over  all  q,  then  the  result  will  be  the  infinite  sequence  of  integro- 
different ial  equations  given  as 

^[ST+T  *Z-l(z'u)  +  22++“l  q>2+l^Z' u)]  +  Et(u)<MZ’u) 

<=/ <iu'ZeS(u' )92(Z,u' )P2(u-u» )  +  80^  ;  f  du' )<p0(z> u*  )h(u'  -  u) 

+  60  ,  /  du'v(u')Sf(u')q>0 (^u');  (2  -  0,  1,  2,  .  .  .,  »)  (4) 

’  *-eff  •'O 

where  8q  2  is  the  Kronecker  delta  defined  here  by  the  properties 


»0,  l  -  1*0 

=  0,  2/0 


(5) 


* 


[ 


Equation  (4)  represents  the  spherical  harmonics  form  of  the  transport 
equation  in  slab  geometry. 
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Slab  Geometry  Spherical  Harmonics  Form  of  Angular  Slowing- 

Down  Density  Due  to  Elastic  Scattering 

Another  important  characteristic  of  slowing  down  neutrons  is  the 
slowing-down  density  Q(r,  u,  j^).  The  slowing-down  density  due  to  elastic 
scattering  Q(r,  u,  If)  is  defined  as  the  number  of  neutrons  which  pass  the 
lethargy  mark  u  per  unit  time  per  unit  phase  volume  (dr  dfi).  This 
quantity  is  easily  derivable  as 

_  DO 

Q(r,u,n)  =  f  du'  f  du"  f  an 'S  (u*  )cp(r,ujfi'  )P(u"  -  u/(iL)  (6) 
•'O  */u 

If  the  slowing-down  density  varies  only  along  the  Z-direction  and 
is  azimuthally  symmetric  about  that  axis,  if  equation  (6)  is  integrated 
over  the  azimuth  angle  <p,  if  the  expressions  for  cp(Z,  u,  p)  and 
P(u-u,'|1l),  given  by  equations  (3a)  and  (3b),  respectively,  are  sub¬ 
stituted  into  (6)  and  then  integrated  over  di2',  and  if  the  resulting 
equation  is  multiplied  by  a  Legendre  polynomial  PK^  and  integrated 

over  all  u,  then  the  result  will  be  the  moments  Q^(Z, u)  of  the 
slowing-down  density  Q(Z,  u,j?);  that  is, 

/*U  /•" 

Qz(Z,u)  -  f  d.f2Q(Z,u,iT)Pz(n)  =  J  du  J  du"Ses(u')q>z(Z,u')Pz(u"  -  u' ) 

( l  =  0,  1,  2,  .  .  .,  co)  (7) 

Equation  (?)  represents  the  spherical  harmonics  form  of  the  slowing-down 
density  in  slab  geometry. 


Lethargy- Dependent  Consistent  PI  Equations  in  Slab  Geometry 

To  obtain  the  PI  equations  from  equation  (4)  the  assumption  is  now 
made  that  (p^Cz, u)  a  0  for  7  ^  2,  which  implies  for  equation  (7)  that 
u)  a  0  for  l  t.  2  also.  The  lethargy-dependent  consistent  PI 
equations  can  now  be  obtained  by  assuming  the  following  spatial  and 


lethargy  dependence  for  9q(Z, u)  and  cp^(z,  u): 

cpQ(Z.,u)  =  cp(u)G[Bm(u)z] 

(8a) 

q^fou)  =  J(u)F[Bm(u)z] 

(8b) 
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In  addition,  the  material  buckling  B^u)  is  assumed  to  be  a  slowly- 
varying  function  of  lethargy.  Therefore,  the  integral  terms  arising  in 
equation  (4)  can  be  evaluated  by  taking 


► 


i 
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cp0(Z,u')  =  cp(u’)G[\(u)z] 

(8c) 

<Pi(z,uO  ^  J(u')?[\(u)z] 

(Qd) 

The  lethargy  and  space  variation  for  the  2  =  0, 1 
follows  from  equations  (8a)  and  (8b)  and  is 

terms  of  equation  (7) 

Qq(Z,u)  =  Q(u)5[\(u)Z] 

(8e) 

Qq(Z,u)  *  P(u)P[^(u)z] 

(8f) 

where  G  and  P  are  spatial  functions  and  Bm(u) 
lethargy-dependent  material  buckling.  In  addition, 
Q  and  P  must  satisfy  the  conditions 

is  defined  as  the 
the  spatial  functions 

^(o[%(u)z])  -  -  B„(-i)F[rB(a;z] 

(9a) 

(9b) 

With  the  assumptions  implied  by  equations  (8a)  to  (8d)  and  (9a)  and  (9b), 
the  spatial  dependence  can  be  divided  out,  and  the  lethargy-dependent 
consistent  PI  equations  obtained  from  equation  (4)  are 

B^u^u)  +  2t(u)cp(u)  =  J  du'Eggtu'Mu'JPofa^') 

+ y  du'Z^n(u'  )cp(u'  )h(u'  -*  u) 

+  f  du'v(u'  JZ-pfu1  )<p(u' )  (10a) 

Aeff  Jq 


together  with 


"  |  +  Et(u)J(u)  =  f  du'Eea(u,)Jr(u')l>i(u-u' )  (i0b) 

h 

Q(u)  =  y*  du '  J  du"Eeg(u'  )cp(u'  )Po(u"  -u')  (lla) 
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P(u) 


au,,Zes(u')j(u,)P-L(u1' 


u1) 


(m) 


In  reference  12  Simon  discusses  some  analytical  solutions  of  the 
lethargy-dependent  consistent  PI  equations  for  hydrogen.  To  obtain  the 
analytical  results,  special  cross-section  variations  were  assumed  for 
hydrogen. 


Elastic  Scattering  Kernel  for  Linearly  Anisotropic 


Center  of  Mass  Elastic  Scattering 


Since  most  light  elements  exhibit  anisotropic  elastic  scattering 
in  the  center  of  mass  coordinate  system  for  high-energy  neutrons,  it  is 
important  uu  include  this  effect  if  accurate  values  of  the  neutron  age 
are  to  be  computed.  Therefore,  by  following  Marshak  (ref.  13)  it  is 
easy  to  show  that  for  linearly  anisotropic  scattering  in  the  center  of 
mass  coordinate  system  the  elastic  scattering  kernel  F(u-u,Vt,)  Is 
given  by 


P(u~u,'nL)du  duL  dcpL 

_  (A  +  I)2  ~ (u—u ' ) 


1  +  3fn 


]>')[• 

7T  ■  5^L(u-u')  -^±1  e  2 


(A  +  P.)  e-(u-u') 

2A  e 


■|(u-u' )  A  -  1  2^u"u' ^ 


2A 

X  du  d|iL  dq^ 


(12) 


where  fp(u)  is  the  average  cosine  of  the  angle  between  the  initial  and 
final  direction  of  motion  in  the  center  of  mass  coordinate  system  of  the 
elastically  scattered  neutron  of  initial  lethargy  u,  A  is  the  ratio 
of  the  mass  of  the  scattering  nucleus  to  that  of  the  neutron,  8  is  the 
Dirac  delta  function,  and  is  an  appropriate  azimuthal  angle  in  the 

laboratory  coordinate  system.  The  moments  Fj(u-u')  of  the  expansion 
of  P(u-u/hl)  into  a  series  of  Legendre  polynomials,  as  in  equation  (3b), 
are  given  by  the  orthogonality  property  of  these  polynomials  as 

Pz(u-u')  =/d%PzGiL)P(U-u,VL)  (13) 

Substituting  equation  (12)  into  equation  (13)  and  performing  the  indi¬ 
cated  integrations  for  l  =  0, 1  will  give 
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Equations  (15a)  and  (155)  reduce  to  their  isotropic  center  of  mass 
elastic  scattering  values  if  f-^u’ )  h  0. 


Differential  Form  of  Consistent  PI  Equations 

Although  equations  (10a)  and  (10b)  are  the  lethargy-dependent  con¬ 
sistent  PI  equations,  a  differential  form  of  these  equations  would  be 
more  suitable  for  computation.  To  obtain  these  differential  equations 
differentiate  equations  (lla)  and  (lib)  with  respect  to  lethargy  u: 

Mil  =  -  du'Ees(u'  )cp(u'  )PQ(u-u' )  +  £es(u)cp(u)  f  du"P0(u"  -  u) 


=  -  /*U  du'Ees(u'  )<p(u’  )P0(u-u' )  +  Ses(u)cp(u)  (16a) 
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SM  „  -f  du'Ees(u')j(u,iP1(u-u')  +  Sea(«)j(u )jf  du^fu"  -  a) 

'  ’A  'Seg(u' )J(u' )P1(u-u’ )  +  n0(u)Ees(u)j(u)  (16b) 

/“ 

du"Pn(u"  -  u)  is  evaluated,  using  equation  (15a)  and 
u  u  /*°° 

is  equal  to  unity.  The  integral  J  du"Pq(u"  -  u)  is  evaluated  using 

*u 

(15b)  and  is  equal  to  hq(u)  -  the  average  cosine  of  the  angle  between 

the  initial  and  final  direction  of  motion,  as  measured  in  the  laboratory 
coordinate  system,  of  the  elastically  scattered  neutrons.  The  first 
integral  on  the  right  side  of  equation  (10a)  can  be  evaluated  using  equa¬ 
tion  (16a),  and  the  integral  term  on  the  right  side  of  equation  (10b) 
can  be  evaluated  using  equation  (16b).  Equations  (10a)  and  (10b)  now 
become 

EmMJ(u)  +  [xa(u)  +  %»]?(>>)  + 


_  f  (u)  r  du'V  (u1  )Ef  (u1  )tf>(u'  )  +  f  du'E.-fu'  )cp(u'  )h(u'  -*  u) 
Keff  Sq  J 

(17a) 

I  «*(.)  +  [St(u>  -  J(u)  +  .  0  (17b) 


where  £a(u)  is  the  macroscopic  absorption  cross  section.  Equations 

(17a)  and  (l7b)  represent  the  basic  lethargy-dependent  equations  of  the 
consistent  PI  approximation  to  the  Boltzmann  transport  equation  for 
slowing -down  neutrons. 


Two  more  equations  are  needed  to  complete  the  PI  description  of 
the  slowing-down  process  as  given  by  equations  (17a)  and  (17b).  These 
are  coupling  equations  between  cp(u)  and  Q(u),  and  J(u)  and  P(u), 
which  can  he  obtained  by  writing  approximate  forms  for  equations  (ll) 
and  (16).  These  additional  equations  are  obtained  in  the  following 
manner.  Assume  that  the  functions  £eg(u' )cp(u' )  and  Ees(u')J(u')  are 

slowly  varying  functions  and  that  they  can  be  approximated  by  the  first 
two  terms  of  a  Taylor's  series  expansion  about  lethargy  u.  That  is, 


Ees(u’ Mu')  =  2es(u)(:p(u)  +  (u’ 

-  u)  cE  [Ees(uMu)] 

(18a) 

Ees(u')j(u')  =  Zes(u)j(u)  +  (u' 

“  u)  H  C2es(u)J(u)] 

(18b) 
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Substitute  equations  (l8a)  into  (lla)  and  (16a)  and  equation  (18b)  into 
(lib)  and  (16b).  Then,  evaluating  the  various  integrals  over  lethargy 
of  Pq(u-u')  and  Pp(u-u')  and  eliminating  the  terms  involving 

^  [£es(u)cP(u)]  and  ^  [2es(u)J(u)]  result  in 

r(u)  +  Q(u)  =  |1(u)Ses(u)cp(u)  (19a) 

p(u)  ^  +  P(u)  =  ^(UJZ^OOJOO  (19b) 

The  parameters  appearing  in  equations  (19a)  and  (19b)  are  defined,  in 
terms  of  the  quantity  t  =  u-u',  as  averages  over  the  distribution 
function  Po(t).  That  is, 


ip  -  t 

(20a) 

2hr  s  t2 

(20b) 

1  5  ) 

(20c) 

2 P  5  t2|iL(t) 

(20d) 

where  the  bar  denotes  an  average  value.  Of  course,  the  various  param¬ 
eters  defined  by  equation  (20)  depend  on  lethargy  since  P0(u-u' )  and 
Pp(u-u')  depend  upon  fq(u). 

Thus,  equations  (17a),  (1Tb),  (l9a),  and  (19b)  form  the  four 
coupled  first-order  ordinary  differential  equations  describing  the 
neutron  slowing -down  process  in  the  consistent  PI  approximation  to  the 
Boltzmann  transport  equation.  Included  are  the  effects  of  linearly 
anisotropic  center  of  mass  elastic  scattering,  inelastic  neutron  scat¬ 
tering,  and  lethargy-dependent  geometric  buckling.  Although  equations 
(17)  and  (19)  apply  for  only  one  nuclear  specie,  the  extension  of  the 
equations  to  include  additional  nuclear  species  follows  directly.  In 
particular  a  set  of  equations  (19)  is  written  for  each  element  with 
appropriate  modifications  of  equations  (17). 


Comparison  of  Lethargy-Dependent  Consistent  PI  Equations 

i 

| 

With  Various  Approximate  Theories 

Equations  (17)  and  (19)  are  fairly  general,  and  it  is  worth  noting 
that  under  various  restrictions  simpler  sets  of  equations,  which  have 
been  used  frequently  in  reactor  analysis,  can  be  easily  obtained.  This 
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will  be  done  for  a  few  typical  cases.  The  first  of  these  is  the  age- 
diffusion  equation  with  no  inelastic  scattering  which  can  be  obtained 
from  (17)  and  (19)  by  setting  2in(u)  s  0  in  (l7a),  dP(u)/du  s  0  in 

(17b).,  dQ(u)/du  =0  in  (19a),  and  dP(u)/du  =  P(u)  =0  in  equation 
(19b). 

The  familiar  Fermi  age-diffusion  equations  result  and  are  given  as 
(ref.  14) 

Bm(u)j(u)  +  Ea(u)cp(u)  +  jf  du'v  (u1  )£f (u1  )cp(u' )  (21a) 


J(u)  =  D(u)Bm(u)cp(u) 

Q(u)  =  l1(u)Ees(u)cp(u) 

where  the  diffusion  coefficient  D(u)  has  been  defined  as 

3{£a(u)  +  ^es(u^  L1  "  HoC«>]} 


(21b) 

(21c) 


(22) 


The  Selengut-Goertzel  approximation  is  basically  the  age-diffusion  set 
given  by  equations  (21a)  to  (21c)  but  with  the  element  hydrogen  treated 
by  equation  (19a)  -  other  elements  still  being  treated  by  equation  (21c). 
To  see  this,,  recall  that  hydrogen  exhibits  isotropic  center  of  mass 
elastic  scattering  collisions  over  energies  of  interest  in  reactors. 

Thus,  £i(u)  and  r(u)  for  hydrogen  are  both  equal  to  unity  for  all 
lethargies  giving 


+  qH(u)  =  SHg(u)cp(u)  (23) 

which  is  the  Selengut-Goertzel  equation  relating  the  hydrogen  slowing- 
down  density  QH(u)  to  the  lethargy-dependent  flux  cp(u)  (p*  348,  ref. 
15).  The  Greuling-Goertzel  approximation  (p.  390,  ref.  15)  includes 
both  equations  (21a)  and  (21b)  of  the  age-diffusion  formalism  but  re¬ 
tains  equation  (19a)  for  all  elements.  Other  approximations  are  also 
possible  but  will  not  be  discussed  because  machine  programs  are  just  as 
easily  formulated  for  the  basic  set  of  equations  given  by  (17)  and  (19) 
as  for  even  the  simple  age-diffusion  set  of  equations. 


Equations  for  Thermal  Group 

The  slowing-down  equations  given  by  (17a),  (17b),  (19a),  and  (19b) 
are  assumed  to  apply  to  some  cutoff  lethargy  Up  corresponding  to  an 

energy  Ep.  The  energy  Ep  is  usually  taken  as  about  5kT  where  T 

is  the  temperature  of  the  moderator  in  °K  and  k  is  the  Boltzmann  con¬ 
stant.  In  the  energy  range  from  0  to  fij  the  neutrons  have  energies 
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comparable  to  the  energies  of  the  thermally  agitated  nuclei  and  can 
gain  as  well  as  lose  energy  upon  making  a  collision.  However,  if  the 
absorption  and  leakage  of  neutrons  are  small,  the  neutron  energy  spec¬ 
trum  will  approach  a  Maxwell-Boltzmann  distribution  in  speed  with  a 
temperature  corresponding  to  that  of  the  moderator.  Under  these  con¬ 
ditions  and  using  the  principle  of  detailed  balance,  the  thermal  equa¬ 
tion  in  the  PI  approximation  will  be,  after  separating  the  spatial 
function  as  before, 


du  /  du'S^Cu'  )cp(u'  )H(u’  -»  u)  (24a) 

j  •S-co 


3  +  2trTJT  “  p^ut) 


(24b) 


where  S, 


laT 


and  E. 


are  the  appropriate  averaged  values  of  the  macro* 


scopic  absorption  and  transport  cross  sections  over  the  thermal  speed 
range.  Thus,  equations  (17a),  (17b),  (19a),  (19b),  (24a),  and  (24b) 
form  the  complete  set  of  consistent  PI  equations  extended  to  include 
many  elements,  which  are  the  basis  for  an  IBM  704  code. 


> 


Discussion  of  Boundary  Conditions,  Fission  Spectrum, 
and  Inelastic  Scattering  Kernel 

If  the  zero  of  lethargy  is  taken  sufficiently  lower  than  the  low¬ 
est  lethargy  in  the  source  spectrum,  the  initial  conditions  for  the 
set  of  consistent  Pi  equations  will  be 

Q(u=0)  =  0,  cp(u=0)  =  0  (25a) 

P(u=0)  =  0,  J(u=0)  =  0  (25b) 


Depending  upon  the  reactor,  a  number  of  possible  source  spectrums  can 
be  used.  For  this  study  the  uranium-235  fission  spectrum,  given  in 
analytic  form  as 


f(u) 


=  0.45270  E  exp 


sinhV2.29  E 


(26) 


has  been  used.  The  quantity  E  is  the  neutron  energy  in  Mev.  The 
spectrum  is  normalized  to  unity;  that  is. 


This  normalization  of  the  source  spectrum  permits  the  evaluation  of  the 
static  criticality  factor  of  the  reactor  as 
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(28) 


The  characteristic  dimension  of  the  reactor  is  needed  to  compute  the 
lethargy-dependent  buckling  B^u).  For  example,  in  the  case  of  a 
spherical  reactor  the  equation  used  for  B^u)  is 

Bm(u)  =  r“+  0.71 

where  r  is  the  radius  and  where 

\r(u)  =  3D(u) 

The  inelastic  scattering  kernel  in  equation  (l7a)  has  been  obtained 
from  the  evaporation  model  of  the  nucleus  as  described  by  Blatt  and 
Weisskopf  (ref.  16).  The  evaporation  model  of  the  nucleus  appears  to 
give  the  gross  features  of  heavy  element  inelastic  scattering.  In  terms 
of  energy,  this  model  gives  for  the  inelastic  kernel  the  following 
result j 

const.  E[e“^E')E]  (31) 

from  the  normalization  requirement  that 

H(E'  -  E)  =  1.0  (32) 

In  equation  (3l)  E'  represents  the  initial  energy  of  the  neutron  and 
E  the  energy  after  an  inelastic  scattering  collision  has  occurred.  The 

quantity  £(Er  )  is  the  reciprocal  nuclear  temperature  in  Mev”-1-  and  is  a 
function  of  the  scattering  nucleus  and  initial  energy  of  the  neutron. 

It  has  been  approximated  by 

h*’ )  -  (4f)1/2  <33> 

with  the  parameter  a  being  experimentally  determined  for  each  nucleus 
of  interest. 


H(E’  -  e)  = 

where  the  constant  is  determined 


/E' 


dE 


(29) 

(30) 


Discussion  of  Neutron  Cross  Sections 

The  lethargy-dependent  consistent  PI  equations  described  along 
with  the  initial  conditions,  source  spectrum,  and  approximate  inelastic 
scattering  kernel  are  integrated  for  a  homogeneous  unreflected  reactor 
system  over  the  microscopic  cross-section  data  of  the  elements  involved 
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in  a  stepwise  fashion  from  107  ev  to  0.025  ev  (or  from  lethargy  u  =  0 
to  u  =  19.8).  This  lethargy  range  is  spanned  by  473  mesh  points  such 
that  the  cross  sections  of  the  various  elements  are  well  represented 
in  the  resonance  energy  range.  The  microscopic  cross-section  data  of 
29  elements  are  represented  in  such  a  fashion  for  use  with  the  IBM  704 
code.  These  cross-section  data  include  aa,  aes,  v<7f,  fp,  and 

where  aa  is  the  microscopic  absorption  cross  section,  aes  is  the 
microscopic  elastic  scattering  cross  section,  vo^  is  the  microscopic 
production  cross  section,  f^  is  the  average  cosine  of  the  angle  between 

the  initial  and  final  direction  of  motion  of  an  elastically  scattered 
neutron  as  measured  in  the  center  of  mass  coordinate  system,  and  a^n 

is  the  microscopic  inelastic  scattering  cross  section.  The  quantity  v 
is  the  average  number  of  neutrons  produced  in  fission,  and  Of  is  the 

microscopic  fission  cross  section.  The  quantity  v  is  taken  to  be 
energy-dependent  according  to  the  data  of  reference  17.  In  the  thermal 
energy  range  the  absorption  and  fission  cross  sections  are  corrected  for 
non-l/v  behavior  as  well  as  for  hardening  of  the  thermal  spectrum,  due 
to  absorption,  using  the  data  of  reference  18.  The  transport  cross  sec¬ 
tion  at  thermal  energy  is  similarly  corrected. 


THEORETICAL  RESULTS  AM)  COMPARISON  WITH  EXPERIMENT 
Age  of  Fission  Neutrons  in  Water 

One  of  the  important  parameters  describing  moderating  materials  is 
the  fission  neutron  age  t.  For  hydrogenous  media  the  age  becomes  espe¬ 
cially  important  because  of  the  well-known  discrepancy  for  water  between 
the  experimental  measurements  and  various  theoretical  calculations.  For 
any  homogeneous  medium  the  lethargy-dependent  consistent  PI  equations 
can  be  solved  for  a  given  set  of  microscopic  cross-section  data,  the 
uranium-235  fission  spectrum,  and  for  a  constant  buckling  B,,].  From 
the  neutron  fluxes  and  currents  generated  as  a  function  of  lethargy,  it 
is  possible  to  compute  the  weighted  fission  age. 

In  this  work,  the  weighted  fission  age  for  water  to  indium  reso¬ 
nance  energy  is  computed  to  be  26.5  square  centimeters.  Oxygen  was  con¬ 
sidered  to  have  linearly  anisotropic  center  of  mass  elastic  scattering. 
For  isotropic  center  of  mass  scattering  for  oxygen,  the  age  is  calcu¬ 
lated  to  be  24.7  square  centimeters.  This  demonstrates  the  importance 
of  precise  calculation  of  the  various  competing  effects  on  the  neutron 
spectrum.  For  comparison  of  the  age  of  26.5  square  centimeters  calcu¬ 
lated  in  this  work,  the  weighted  mean  of  a  number  of  different  types  of 
calculations  is  given  in  reference  19  as  26. Q±0- 5  square  centimeters. 

It  is  also  of  interest  to  note  that  a  recent  measurement  of  the  fission 
neutron  age  in  water  is  27.86  square  centimeters  (ref.  20). 
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Fourier  Transform  of  Slowing -Down  Kernel  in  Water 

Another  important  property  of  homogeneous  moderating  media  can  be 
described  in  terms  of  the  Fourier  transform  K(Btn)  of  the  slowlng-down 

kernel.  Physically,  the  Fourier  transform  K(Bm)  is  the  nonleakage 

probability  of  the  source  neutrons  to  a  given  energy.  The  transform 
K(Bm)  can  be  derived  from  a  solution  of  the  lethargy-dependent  consistent 
PI  equations  for  a  given  source  spectrum  and  a  constant  value  of  the 
buckling  B^. 

Figure  1  is  a  plot  of  the  transform  K(Bm)  as  a  function  of 
for  fission  neutrons  in  water  to  indium  resonance  energy.  Shown  on  the 
curve  are  values  of  the  transform  computed  by  Gwin,  Trubey,  and  Weinberg 
(fig.  1,  ref.  6)  from  the  experimental  slowing-down  kernel  as  measured 
by  Hill,  Roberts,  and  Fitch  (ref.  2l).  The  agreement  over  a  wide  range 
of  buckling  is  seen  to  be  good.  The  value  of  the  fission  neutron  age  in 
water  reported  in  reference  21  is  higher  than  the  recent  value  of  ref¬ 
erence  20  by  about  3  square  centimeters.  This  is  believed  due  to  dif¬ 
ferences  between  the  two  experiments  of  measurements  near  the  source. 

This  difference  near  the  source  of  the  reference  21  measurement  of  the 
slowing-down  kernel  in  water  is  expected  to  have  a  small  effect  on  the 
calculation  of  the  transform  K(Bm). 

Table  I  compares  the  transform  K(Bm)  to  indium  resonance  energy 
for  water  and  two  concentrations  of  uranyl  fluoride  salt  in  water  solu¬ 
tion.  The  differences  as  compared  to  water  are  slight  except  for  extreme 
values  of  the  buckling. 


Calculation  of  Critical  Dimensions  of  Unreflected 
Homogeneous  UOgFg-HgO  Assemblies 

The  lethargy-dependent  equations  of  the  consistent  PI  slowing-down 
model  described  in  the  section  on  theory  will  be  used  to  calculate  the 
parameters  of  critical  assemblies  of  unreflected  uranyl  fluoride  salt 
UO2F2  dissolved  in  water.  Comparison  of  theoretical  results  with  avail¬ 
able  experimental  data  will  also  be  made. 

Before  any  calculations  can  be  performed,  the  physical  properties 
of  uranyl  fluoride  -  water  solutions  must  be  known.  Figure  2  is  a  curve 
of  density  of  uranyl  fluoride  -  water  solutions  as  a  function  of  hydrogen 
atom  concentration  for  a  temperature  of  25°  C  obtained  from  available 
experimental  data  (ref.  22).  The  point  on  the  curve  for  pure  water  is 
for  a  density  of  0.997  gram  per  cubic  centimeter  and  a  hydrogen  atom 
concentration  of  6.66X10^'  atoms  per  cubic  centimeter.  Since  the  uranium 
in  the  solution  is  considered  to  contain  93.2  percent  by  weight  of  the 
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uranium-235  isotope,  the  data  of  figure  2  can  be  used  to  derive  a  curve, 
figure  3,  of  the  uranium-235  atom  concentrations  IT28  j_n  atoms  per 
cubic  centimeter  as  a  function  of  the  hydrogen  atom  concentration 
in  atoms  per  cubic  centimeter.  On  figure  3  are  also  indicated  various 
values  of  R  defined  to  be  the  ratio  of  the  hydrogen  atom  concentration 
to  uranium-235  atom  concentration,  R  =  N^/N^.  Various  important  char¬ 
acteristics  of  the  W2'S'2~^2(~>  critical  assemblies  will  be  given  in  terms 
of  this  parameter  R. 

Thus,  specification  of  the  solution  density  for  the  UOgFg-HgO  solu¬ 
tion  reactor  in  question  gives  the  hydrogen  atom  concentration  from  the 
curve  given  by  figure  2.  Then,  from  figure  3  the  uranium-235  concen¬ 
tration  can  be  determined.  The  atom  densities  for  the  other  constituents 
of  the  solution  are  determined  from  the  following  relations: 


Uranium-238:  N28  =  0.0730  N25 

(34a) 

Fluorine:  =  2(W25  +  N28) 

(34b) 

Oxygen:  N^2  «=  i 

(34c) 

By  using  the  atom  concentration  for  the  various  constitutents  of 
the  UC^Fg-HgO  solutions  in  the  manner  described,  calculations  were  per¬ 
formed,  using  the  lethargy-dependent  consistent  PI  equations  over  the 
473  lethargy  points  of  the  microscopic  cross-section  data  extending  from 

10 /  ev  to  0.025  ev,  over  a  wide  range  of  solution  concentrations.  In 
terns  of  the  parameter  R,  the  solution  concentrations  were  varied  from 
R  =  25  to  700.  As  a  result  of  these  many  calculations,  the  static 
criticality  factor  Is  obtained  for  a  given  hydrogen  to  uranium- 

235  atom  concentration  (R)  as  a  function  of  the  buckling  of  the 

system.  The  buckling,  although  capable  of  being  varied  with  the  lethargy 
in  the  IBM  704  code,  was  considered  to  be  lethargy-independent  for  these 
calculations. 

As  examples,  the  lethargy-dependent  flux  cp(u)  is  plotted  as  a 
function  of  lethargy  u  for  critical  reactors  having  solution  con¬ 
centration  values  of  R  =  40  and  472.  The  source  term  for  the  cal¬ 
culations  corresponds  to  1  fission  neutron.  The  flux  for  R  =  40,  a 
system  heavily  loaded  with  uranium-235,  is  shown  in  figure  4.  In  the 
fission  spectrum  lethargy  range  from  about  lethargy  u  =  0  to  u  =  6 
a  resonance  structure  in  the  flux  is  noted  which  corresponds  to  various 
high-energy  scattering  resonances  of  oxygen.  In  the  resonance  energy 
range  from  lethargy  u  =  6  to  u  =  18.2,  the  decrease  in  the  flux  with 
lethargy  occurs  because  of  the  absorption  by  the  uranium-235  nuclei. 

The  dips  in  the  flux  correspond  to  absorption  resonances  in  the  uranium- 
235.  The  effect  on  the  flux  of  the  approximately  l/v  capture  is 
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evident  from  lethargy  u  =  15.0  to  18.2.  The  flux  for  R  =  472,  a 
system  lightly  loaded  with  uranium-235,  is  shown  in  figure  5.  In  the 
lethargy  range  from  u  =  0  to  6  the  flux  for  R  =  472  appears  roughly 
the  same  as  for  the  R  =  40  system.  In  the  lethargy  range  from  u  =  6 
to  u  =  18.2  the  flux  is  constant  from  about  lethargy  u  =  7  to 
u  =  12.0.  From  lethargy  u  =  12.0  to  18.2  the  effects  of  uranium-235 
absorption  resonances  and  l/v  capture  are  evident  but  are  not  as 
(drastic  as  for  the  heavily  loaded  R  =  40  reactor  system. 

Figure  6  is  a  plot  of  the  static  criticality  factor  Ke-fj  for  un¬ 
reflected  UO2F2-H2O  assemblies  as  a  function  of  the  parametric  R  and 
the  buckling  B^.  Curves  of  the  criticality  factor  varying  from  0.98 
to  1.12  are  shown.  All  the  curves  exhibit  a  broad  minimum  in  the  buck¬ 
ling  over  values  of  R  from  40  to  70.  For  the  static  criticality 
factor  equal  to  unity  the  radius  of  the  equivalent  unreflected  sphere, 
which  includes  the  extrapolation  distance,  varies  from  about  18. 25  to 
25  centimeters  over  the  range  of  the  parameter  R. 


Comparison  of  Calculated  Critical  Dimensions  with  Experimental 
Values  for  Unreflected  UQ2F2-H2O  Assemblies 

In  order  to  compare  calculated  critical  dimensions  with  experi¬ 
mentally  determined  values,  the  effect  of  the  aluminum  containment 
vessels  on  measured  critical  dimensions  of  unreflected  UO2F2-H2O  assem¬ 
blies  must  be  known.  The  worth  of  the  aluminum  bottom  plate  for  the 
NASA  experiments  was  measured  in  terms  of  the  change  in  critical  height 
as  a  function  of  the  thickness  of  the  bottom  plate  for  a  12-inch- 
diameter  cylindrical  tank.  The  solution  concentration  R  was  475. 

These  results  are  shown  in  figure  7  where  it  is  seen  that,  for  example, 
a  1-centimeter-thick  bottom  plate  is  worth  about  0. 3  centimeter  of  solu¬ 
tion  height  for  this  particular  reactor.  It  is  expected  that  the  bottom 
plate  worths  for  other  values  of  R  would  be  nearly  the  same  as  for  the 
values  for  R  =  475  since  the  fraction  and  spectrum  of  neutrons  leaking 
from  these  do  not  differ  greatly.  It  is  estimated  that  the  reported 
critical  height  for  the  NASA  experiments  could  be  too  low  by  a  maximum 
of  about  0.5  centimeter  since  the  wall  thickness  of  the  aluminum  cylinder 
is  l/8  inch.  Thus,  because  of  the  expected  smallness  of  this  correction, 
none  of  the  Oak  Ridge  and  NASA  criticality  data  have  been  corrected  for 
the  criticality  effect  of  the  aluminum  containment  vessels. 

Since  figure  G  is  based  on  an  energy- independent  buckling,  this 
implies  an  energy- independent  extrapolation  distance  5.  This  value  of 
5  can  be  deduced  from  experimental  criticality  data  on  unreflected 
cylinders  if  the  top,  bottom,  and  side  extrapolation  distances  are 
assumed  to  be  equal.  Then,  comparison  with  a  wide  range  of  reactor 
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dimensions  and  solution  concentrations  will  establish  the  validity  of  a 
constant  value  for  5.  Figure  8  shows  theoretical  curves  based  on  values 
of  S  =  2.5,  3.0,  and  3.5  centimeters  for  unreflected  12-inch-diameter 
cylindrical  reactors  of  UOgFg-HgO  solution  in  which  the  critical  height 

is  plotted  as  a  function  of  the  parameter  R.  These  curves  are  based  on 
the  following  equation  for  the  critical  height  h  of  the  cylinders: 


where 


(35a) 


(35b) 


and  r  is  the  physical  radius  of  the  cylinders.  For  a  given  value  of 
R,  is  evaluated  from  figure  6  for  a  Kg.^  of  unity,  and  then,  since 

the  values  of  r  and  S  have  been  fixed,  the  physical  height  h  is 
given  by  equations  (35a)  and  (35b).  It  is  evident  from  figure  8  that 
these  critical  assemblies  are  quite  sensitive  to  the  value  of  the  ex¬ 
trapolation  distance  which  is  used.  For  the  12 -inch-diameter  cylindrical 
reactors  considered  the  value  of  5  =  3.0  centimeters  seems  to  predict 
the  experimental  critical  heights  listed  in  table  II  quite  well  as  a 
function  of  the  atom  ratio  parameter  R.  This  insensitivity  of  5  with 
R  is  probably  due  to  the  fact  that  the  hydrogen  atom  density,  and  hence 
the  leakage  fraction,  changes  but  little  over  the  range  of  R  considered. 
Therefore,  it  was  decided  to  use  this  value  of  6  =  3. 0  centimeters  for 
the  extrapolation  distance  for  all  values  of  the  parameter  R  for  the 
various  diameters  of  the  unreflected  cylindrical  critical  assemblies  of 
UOgFg-HgO.  Reference  6  reports  using  values  of  8  =  2.5  to  3.7  for  the 
empirical  kernel  method  of  calculation. 

The  curves  of  figure  9  represent  calculations  of  the  critical 
heights  of  cylindrical  reactors  for  various  diameters  and  as  a  function 
of  the  parameter  R.  These  calculations  are  based  on  the  theoretical 
buckling  curve  given  as  figure  6  and  the  assumed  constant  value  of  the 
extrapolation  distance  8  =  3.0  centimeters.  The  experimental  points, 
listed  in  table  II,  represent  data  from  references  2,  3,  4,  and  the 
NASA  solution  reactor.  The  agreement  for  various  diameter  reactors  as 
a  function  of  R  is  quite  good  between  the  theoretical  model  described 
in  this  work  and  the  various  experiments.  Thus,  this  provides  the 
verification  of  the  general  buckling  curve  (fig.  6)  and  of  the  constant 
value  of  3.0  centimeters  for  the  extrapolation  distance  for  these  un¬ 
reflected  cylindrical  reactors. 
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Calculation  of  Nonthermal  Fission  Fractions,  Leakage  Fractions, 
Neutron  Age,  and  Two-  and  Three-Group  Macroscopic  Cross 
Sections  as  a  Function  of  Solution  Concentration 


for  UOgFg-HgO  Assemblies 

In  solving  for  the  critical  dimensions  of  unreflected  UOgFg-HgO 

reactors  using  the  consistent  PI  equations,  the  lethargy-dependent  flux 
cp(u)  and  current  J(u)  are  generated.  Since  all  the  information  about 
a  given  reactor  is  contained  by  the  flux  and  current,  it  is  possible  to 
derive  many  useful  quantities  pertaining  to  the  UOgFg-HgO  reactors. 

Figures  10,  11,  and  12  present  some  of  this  information. 


In  figure  10  is  plotted  the  nonthermal  fission  fraction  as  a  func¬ 
tion  of  the  parameter  R  for  various  cutoff  energies.  For  a  heavily 
loaded  UO^Fg-HgO  system  having  R  =  25,  about  65  percent  of  the  non¬ 


thermal  fissions  occur  between  energies  of  10  ev  to  0.125  ev,  while 
about  43  percent  of  the  total  fissions  occur  between  energies  of  10^  ev 
to  1.44  ev.  For  a  UO2F2-H0O  reactor  having  R  =  700  the  corresponding 


numbers  are  about  6  and  3  percent. 


Figure  11  is  a  plot  of  the  fraction  of  the  neutrons  which  escape 
from  critical  UOgFg-^O  systems  as  a  function  of  R.  This  fraction  is 

nearly  constant  over  the  range  of  R  from  25  to  about  300,  varying  be¬ 
tween  44  and  47  percent.  For  larger  values  of  R  this  function  de¬ 
creases  and  is  about  33  percent  at  R  =  700. 

Figure  12  is  a  plot  of  the  fission  neutron  age  r  as  a  function 
of  the  parameter  R  both  for  the  indium  resonance  detector  energy  of 
1.44  ev  and  the  thermal  cutoff  energy  of  0.125  ev.  From  about  R  =  100 
to  R  =  700  the  age  is  nearly  constant.  For  values  of  R  less  than 
100  the  age  increases  because  of  the  decreased  hydrogen  atom  concen¬ 
tration.  In  the  limit  as  R  -*  »  the  1.44-ev  curve  would  approach  the 
fission  neutron  age  in  water. 


The  neutron  fluxes  generated  in  solving  for  the  critical  dimensions 
of  unreflected  UO2F2-H2O  systems  can  also  be  used  as  a  weighting  function 

for  obtaining  two-  and  three-group  macroscopic  cross-section  parameters. 
These  two-  and  three-group  cross  sections  can  be  used  in  calculating 
reflected  UO2F2-H2O  solution  reactors. 
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Both  three-  and  two-group  constants  were  generated  as  a  function  of 
R.  Group  I  of  the  three-group  split  extends  from  u  =  0  to  u  =  6.1, 
group  II  from  u  =  6.1  to  u  =  18.2,  and  the  thermal  group  is  centered 
about  the  lethargy  of  19.8  corresponding  to  the  room  temperature  kT  of 
0.025  ev.  For  the  two-group  split,  the  thermal  group  is  the  same  as 
that  of  the  three-group  split,  and  the  fast  group  (group  l)  extends  from 
lethargy  u  =  0  to  u  =  18.2.  For  the  few-group  calculations  performed 
with  a  one-dimensional  multigroup  IBM  704  code,  the  parameters  needed 
are  the  group  diffusion  coefficient  DG,  group  macroscopic  absorption 
cross  section  group  macroscopic  production  cross  section  vE^,  and 

the  group  removal  cross  section  These  group  parameters  are  defined 

over  the  fluxes  and  currents  generated  for  a  given  value  of 


R  as 


/ 


j(u)du 


"group 


/ 
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(36a) 


Bni(u)cp(u)du 


£x  = 


"group 


sa(u)tp(u)du 
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(36b) 


cp  ( u )  du 


vE® 


/ 


group 


v  (u)Ef(u)cp(u)du 


"group 


/ 


cp(u)du 


(36c) 


"group 


r« 

The  quantity  E^  is  evaluated  by  making  a  neutron  balance  over  the 

group  in  question.  Over  the  range  of  buckling  considered,  the  group 
constants  were  not  sensitive  to  the  choice  of  buckling  Bpj. 


The  few-group  constants  obtained  by  the  method  described  are  shown 
1r.  figures  13,  14,  15,  and  16  for  both  the  three-  and  two-group  splits 
over  the  range  of  R  values.  The  diffusion  coefficient  D  is  plotted 
in  figure  13  and  is  relatively  constant  over  the  entire  range  of  R  for 
most  of  the  fast  groups;  however,  the  effect  of  uranium-235  absorption 
is  seen  for  all  the  R  values  for  the  thermal  neutron  group.  The  ab¬ 
sorption  cross  section  E^  is  shown  in  figure  14  and  varies  strongly 
with  the  atom  ratio  R.  The  production  cross  section  vE^  is  shown  in 

figure  15  and  exhibits  similar  behavior  as  the  absorption  cross  section. 
The  removal  cross  section  £q  is  shown  in  figure  16  and  increases  as 
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the  parameter  R  increases.  Table  III  gives  the  corresponding  three- 
and  two-group  data  for  water  at  room  temperature. 


o 


i 


Comparison  of  Calculated  and  Experimental  Critical  Dimensions 


for  Water-Reflected  UO2F2-H2O  Reactors 


In  order  to  perform  theoretical  calculations  on  fully  reflected 
cylinders,  the  reflector  savings  must  be  established  as  a  function  of 
reflector  thickness  and  concentration  parameter  R.  Using  the  three- 
group  parameters  plotced  in  figures  13  to  16  for  various  values  of  R 
and  the  water  constants  listed  in  table  III,  the  critical  core  radius  r 
was  calculated  for  spherical  reactors  of  UOgFg-I^O  solutions  with  various 
thicknesses  of  water  using  an  IBM  704  diffusion  code.  The  reflector 
savings  RS  is  defined  as 


RS 


r 


(37) 


2 

where  for  a  given  value  of  R  the  value  of  is  obtained  from  fig¬ 

ure  6.  The  results  of  these  calculations  are  shown  in  figure  17.  The 
curves  show  a  slight  increase  in  the  reflector  savings  with  R.  For  the 
completely  reflected  case  (GO  cm  of  water)  the  value  of  reflector  savings 
is  6.5  centimeters  at  R  =  25  and  S.O  centimeters  at  R  =  700.  For  the 
range  of  R  values  considered  the  reflector  savings  reached  its  maximum 
value  at  about  15  centimeters  thickness  of  water. 


The  critical  heights  h  of  fully  reflected  cylindrical  assemblies 
were  computed  for  cores  of  radius  r  as  a  function  of  the  parameter  R 
using  the  critical,  buckling  cui’ve  of  figure  6  and  the  infinite  water 
reflector  savings  curve  of  figure  17.  These  results  are  shown  in  fig¬ 
ure  18  and  were  calculated  with  the  following  two  formulas: 


(38a) 


(38b) 


Experimentally  obtained  data  listed  in  table  IV  are  also  shown  on  the 
curves.  The  agreement  is  generally  good.  The  dashed  curves  were  com¬ 
puted  by  arbitrarily  increasing  the  reflector  savings  by  0.5  centimeter 
for  all  values  of  R.  The  two  sets  of  curves  bracket  the  experimental 
data  quite  well  and  serve  to  indicate  the  sensitivity  of  the  calculations 
with  small  changes  in  the  value  of  the  reflector  savings. 


1 


i 
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CONCLUSION 

The  lethargy-dependent  equations  of  the  consistent  PI  approxima¬ 
tion  to  the  Boltzmann  transport  equation  for  slowing  down  neutrons,  as 
developed  in  this  present  study,  have  been  used  as  the  basis  of  an  IBM 
704  computer  program.  This  program  has  been  used  to  compute  the  age  of 
fission  neutrons  in  water  as  26.5  square  centimeters  as  compared  to  a 
recent  experimental  value  given  as  27.86  square  centimeters.  The 
Fourier  transform  of  the  slowing-down  kernel  for  water  to  indium  reso¬ 
nance  energy  calculated  here  compares  well  with  the  Fourier  transform  of 
the  kernel  for  water  measured  by  Hill,  Roberts,  and  Fitch.  This  method 
of  calculation  has  also  been  applied  to  uranyl  fluoride  -  water  solution 
critical  assemblies.  Theoretical  results  established  for  both  unre¬ 
flected  and  fully  reflected  critical  assemblies  compare  well  with  avail¬ 
able  experimental  data. 


Lewis  Research  Center 

National  Aeronautics  and  Space  Administration 
Cleveland,  Ohio,  August  22,  1961 
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APPENDIX  -  SYMBOLS 

A  ratio  of  nuclear  mass  to  neutron  mass 

a  inelastic  scattering  parameter  defined  by  eq.  (33) 

Bm(u)  lethargy-dependent  material  buckling  defined  by  eqs.  (8) 

and  (9) 

o  transverse  buckling  for  cylindrical  geometry 

W  D(u)  diffusion  coefficient  at  lethargy  u  defined  by  eq.  (22), 

cm 

D1^  neutron  diffusion  coefficient  for  group  G  defined  by 

eq.  (36a),  cm 

— > 

.  dr  spatial  volume  element 

du  element  of  lethargy 

dv  veJ.ocity  volume  element 

E  neutron  energy 

Eqi  thermal  cutoff  energy  separating  thermal  energy  range 

from  fast  energy  range 

F  spatial  function  defined  by  eq.  (8b) 

f(u)du  fraction  of  source  neutrons  in  lethargy  interval  du 

about  lethargy  u 

f-j(u)  average  cosine  of  angle  between  initial  and  final  direc¬ 

tion  of  motion,  measured  in  the  center  of  mass  coordi¬ 
nate  system,  for  elastically  scattered  neutrons  of 
lethargy  u 

G  spatial  function  defined  by  eq.  (8a) 

H(v'  ->  v)drdv  probability  that  a  neutron  in  phase  volume  element 

(dr  dv' )  will,  after  an  inelastic  scattering  collision, 
be  in  the  phase  volume  (dr  dv) 

h  physical  height  of  cylindrical  reactor 

J(u)  lethargy-dependent  neutron  current  defined  by  eq.  (17b) 
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K(3n) 

Kcff 

k 

N 

n(r, v)dr  dv 

P(u) 

P(v'  -*v)dr  c 
PZ(u-u') 

P  l(») 

Q(r,  u,ft) 

Q(u) 

Q*(Z,u) 

R 

RS 

r 

r 

S(r,v)dr  dv 
T 


Fourier  transform  of  neutron  slowing-down  kernel 
static  criticality  factor  for  a  reactor 
Boltzmann  constant 

concentration  of  an  element,  atoms/cm'5 

probable  number  of  neutrons  in  space  volume  element  dr 
about  position  r  and  in  velocity  element  dv  about 
velocity  v  per  unit  time 

anisotropic  slowing-down  density  defined  by  eq.  (19b) 

probability  that  a  neutron  in  the  phase  volume  element 
(dr  dv' )  will,  after  an  elastic  scattering  collision, 
be  in  the  phase  volume  (dr  dv) 

lethargy-dependent  coefficient  of  index  l ,  in  expansion 
of  elastic  scattering  kernel  into  a  series  of  Legendre 
polynomials,  defined  by  eq.  (13) 

Legendre  polynomial  or  argument  p  and  index  2 

slowing-down  density  defined  as  number  of  elastic  scat¬ 
tered  neutrons  which  pass  lethargy  mark  u  per  unit 
time  per  unit  phase  volume  (dr  dfi) 

isotropic  slowing-down  density  defined  by  eq.  (19a) 

coefficient  of  index  2  in  expansion  of  slowing-down 
density  into  a  series  of  Legendre  polynomials  defined 
by  eq.  (7) 

ratio  of  hydrogen  atom  concentration  to  uranium  atom  con¬ 
centration  for  the  UOgFg-HgO  solution  reactors 

reflector  savings  defined  by  eq.  (37) 

physical  radius  of  cylindrical  reactors 

position  vector 

source  of  neutrons  due  to  fission  which  are  produced  in 
phase  volume  element  (dr  dv) 

moderator  temperature,  °K 
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t  contraction  for  lethargy  difference  u-u' 

u  neutron  lethargy  defined  in  terms  of  neutron  energy  E 

by  equation  u  =  ln(EQ/E),  Eq  usually  taken  as  10  ev 

uj  cutoff  lethargy  corresponding  to  cutoff  energy  Ep 

v  neutron  speed  corresponding  to  energy  E 

■  ^ 

v  neutron  velocity 

Z  coordinate  axis 


r 


8 

80,  l 

5 

\r 


^  L 


Uq 


^1 


v(E) 


h 

p 


quantity  related  to  mean  squared  lethargy  increase  per 
collision,  defined  by  eq.  (20b) 

Dirac  delta  function  used  in  eq.  (12) 

extrapolation  distance  defined  by  eqs.  (35a)  and  (35b) 

Kronecker  delta  defined  by  eq.  (5) 

reciprocal  nuclear  temperature  defined  by  eq.  (33) 

neutron  transport  mean  free  path  defined  by  eq.  (30) 

cosine  of  angle  between  Z-axis  and  direction  of  motion 
of  neutron  in  planar  geometry 

cosine  of  angle  between  initial  and  final  directions  of 
motion  measured  in  laboratory  coordinate  system  of  a 
neutron  undergoing  an  elastic  scattering  collision 

average  cosine,  measured  in  laboratory  coordinate  system, 
for  elastically  scattered  neutrons 

mean  of  the  product  of  lethargy  increase  in  a  collision 
times  angle  of  scattering,  defined  by  eq.  (20c) 

average  number  of  neutrons  produced  per  fission  occurring 
at  energy  E 

average  increase  in  lethargy  per  collision 

quantity  related  to  mean  of  the  product  of  squared 
lethargy  increase  per  collision  times  angle  of  scat¬ 
tering,  defined  by  eq.  (20d) 
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macroscopic  neutron  absorption  cross  section,  cm"'1' 

macroscopic  neutron  elastic  scattering  cross  section,  cm"-*- 

macroscopic  neutron  fission  cross  section,  cm“^ 

macroscopic  neutron  inelastic  scattering  cross  section, 
cm“l 

macroscopic  neutron  total  cross  section,  £a  +  Ees  +  E^n, 
cm--*- 

macroscopic  neutron  absorption  cross  section  for  group  G, 
defined  by  eq.  (36b),  cm"'1 

macroscopic  neutron  fission  cross  section  for  group  G, 
defined  by  eq.  (36c),  cm"'*' 

macroscopic  neutron  removal  cross  section  for  group  G, 
cm"-’- 

microscopic  neutron  absorption  cross  section,  cm^ 

O 

microscopic  neutron  elastic  scattering  cross  section,  cm 

p 

microscopic  neutron  fission  cross  section,  cm 

O 

microscopic  neutron  inelastic  scattering  process,  cm^ 

neutron  age,  defined  as  one-sixth  the  mean  squared 
slowing-down  distance  to  energy  E  of  a  neutron 
originating  from  a  point  source  emitting  neutrons  of 
energy  Ef  in  an  infinite  medium 

neutron  flux,  nv 

unit  vector  in  direction  of  motion  of  neutron 
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TABLE  I.  -  FOURIER  TRANSFORM  K(\)  OF  SLOWENG-DOWN 
KERNEL  FOR  FISSION  NEUTRONS  IN  WATER  AND  IN 
U02F2-H20  SOLUTIONS  WITH  ATOM  RATIOS 
OF  HYDROGEN  TO  URANIUM-235,  R, 


EQUAL  TO  50,  150 


K(Bj 

Bm 

h2o(r  =  ») 

R  =  150 

R  =  50 

0 

1.0 

1.0 

1.0 

. 03142 

.9745 

.9748 

.9752 

.06283 

.9047 

.9057 

.9069 

. 12566 

.6944 

.  6970 

.7000 

.20944 

.  4188 

.  4220 

.4263 

. 26180 

.2881 

.2922 

.2970 

.31416 

,1923 

.1963 

.2013 

.  3927 

.0995 

.1030 

.1079 

.5236 

.0205 

.0310 

.0349 

.  7854 

.0005 

.  0013 

.0028 
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TABLE  III.  -  THREE-  AND  TWO-GROUP  CONSTANTS  FOR  WATER  REFLECTOR 


Lethargy, 

u 

Diffusion 

coefficient, 

dg 

Absorption 
cross  section, 

Removal 
cross  section, 

0-10.2 

1.274 

0.00066 

0.04598 

0-6.1 

1.722 

.00012 

.  07837 

6.1-18.2 

.612 

.00145 

. 11392 

19.8 

.1329 

.01962 

TABLE  IF.  -  EXPERIMENTAL  CRITICALITY  DATA  FOR  U02F2-H20  REACTORS 


FULLY  REFLECTED  BY  AN  INFINITE  WATER  REFLECTOR 


Cylinder 

diameter, 

in. 

Atom  ratio, 

R  «  nH/jjU-235 

Critical 

height, 

cm 

Reference 

8 

26.2 

21.5 

2,  I 

).  69 

29.9 

20.7 

52.9 

19.5 

58.8 

20.5 

99.5 

22.4 

192.0 

2Q.1 

1 

I 

290.0 

40.1 

T 

43.2 

18.6 

4,  I 

>.  72 

10 

52.9 

13.4 

2,  I 

>,  69 

328.7 

22.4 

J 

499 

35.2 

43.2 

12.5 

4,  I 

).  72 

51.5 

11.4 

4,  I 

).  72 

15 

52.9 

7.90 

2,  I 

).  69 

56.7 

8.50 

221.0 

11.30 

499.0 

16.90 

755.0 

27.10 

issssmtsiBmiMBZ  s  a  w—fl 

BBBBBBSS  l!  !■■■■■  I  :-B— 
BMSSMI8SSSSS5HglBB8aHliM 

■  !  1  !■ 

lllllllBliiiiiilSliiliiliililiiiii 


■■■issssasaaBBiBneinsBiiB 
MBHgwBBaaaHBgMBM 

■  ;M  Si  9  Blil 

Eli  i!  SS  ■!  Hi  HS 
■BiiiiBiBiEggaaaaBBBaBEBiaBaBi 

I  ;  35 

■iiiin 

in  mammsssssammmmm 

SB  !llii  IIIH  ^  !! 

H  iiiiisiiiiii  &  a  iiii 
■9ElBSllBSiSl:  ;  il  IE 

■SiHMnBHMHannpaipa 

arcs*  H 

ImmammmmmmammmmmmmmammmmmemmammmmmmmM . . 

■■■■■■■HBMBnBaaapB 
HMMI  «  ■  -  '  -  MHMMMM 

It  ipamnni  ■■■■allii  saunlMniigt 
KmMHHmmmmummamwmmmBmMwSmmmwmamaSBBmBm 

^nmnagMSTHsai  iii 
S:1  iriii  !H 


Uran.1  um-235  atom  concentration  in  UC2?2""2^  solutions.  Uranium  enriched  zc  C.932  cranium-235;  water  at  25°  C;  density,  O.S97Q  gram 
centimeter. 


neutron.  Therxal  neutron  flux  *  O.*704j  /  f(u)du  «  12.582. 
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Atom  ratio,  R  *  N*Vn^"235 


Figure  9.  -  Critical  height  plotted  against  atom  ratio  N^/N^-235  for  various 
diameters  of  unreflected  UO2F2  solution  reactors.  (Theoretical  curves  based 
on  total  buckling  from  figure  6  and  extrapolation  distance  6  *■  3.0  era.) 
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Figure  10.  -  Fraction  of  nonthermal  fissions  to  total  fissions  plotted  against  atom  ratio 
JJH/hU-235  for  several  cutoff  energies  for  UOzFz-EzO  solution  reactors. 
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.44  ev)  and  to  thermal  cutoff  energy  (0-125  ev  =  5  kT)  for  UO2F2-H2O  solution 


Figure  13.  -  Two-  and  three-group  diffusion  coefficients  D  plotted  against  atom  ratio 
for  the  UO2F2  solution  reactors. 
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Figure  14-  -  Two  and  three-group  absorption  cross  sections  plotted  against  atom  ratio 

jjH/jjU-235  for  UO2F2  solution  reactor. 
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